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Abstract. The first galaxies form due to gravitational collapse of primordial halos. During this 
collapse, weak magnetic seed fields get amplified exponentially by the small-scale dynamo - a 
process converting kinetic energy from turbulence into magnetic energy. We use the Kazantsev 
theory, which describes the small-scale dynamo analytically, to study magnetic field amplification 
for different turbulent velocity correlation functions. For incompressible turbulence (Kolmogorov 
turbulence), we find that the growth rate is proportional to the square root of the hydrodynamic 
Reynolds number, Re^''^. In the case of highly compressible turbulence (Burgers turbulence) the 
growth rate increases proportional to Re'^^. With a detailed chemical network we are able 
to follow the chemical evolution and determine the kinetic and magnetic viscosities (due to 
Ohmic and ambipolar diffusion) during the collapse of the halo. This way, we can calculate the 
growth rate of the small-scale dynamo quantitatively and predict the evolution of the small- 
scale magnetic field. As the magnetic energy is transported to larger scales on the local eddy- 
timescale, we obtain an estimate for the magnetic field on the Jeans scale. Even there, we find 
that equipartition with the kinetic energy is reached on small timescales. Dynamically relevant 
field structures can thus be expected already during the formation of the first objects in the 
Universe. 
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1. Introduction 



The present-day Universe is strongly magnetized, which is shown by observations of 
various astrophysical phenomena like jets and outflows caused by magnetic fields. The 
problem however is that theory predicts the first seed field to be extremely weak. Thus, 
we need to explain how these fields are amplified and reach the observed present-day 
large field strengths. 

A process that can amplify seed fields on very short timescales is the small-scale or 
turbulent dynamo. The small-scale dynamo converts kinetic energy from turbulence into 
magnetic energy very efficiently. It strongly depends on the nature of turbulence and here 



we explore the range reaching from incompressible Kolmogorov turbulence ( Kolmogorov 
1941| ) to highly compressible Burgers turbulence ( jBurgers ||1948 ). 



The first point in time where the dynamo could operate is during the formation of the first 
stars at a redshift of about 20. Here the turbulence in the primordial halo is generated by 
accretion ( Greif et al. [2008 Wise et al. [2008 Latif et al. [2012 ) and the weak seed fields 
created in the very early Universe or by battery processes are amplified by the small-scale 
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dynamo (Schleicher et al. 2010 Federrath et al. 11201 la Sur et al. 2012 Peters et al. 



20121. 



The outline of this article is as follows: In Sec.|2]we discuss the basic assumptions of our 
model including the strength of the magnetic seed fields, the description of turbulence 
and the chemical and thermal evolution of primordial gas during collapse. Furthermore, 
we calculate the typical MHD quantities like the hydrodynamic and magnetic Reynolds 
number. Sec. |3]is on the kinematic growth phase of the small-scale dynamo. We discuss 
the influence of different types of turbulence on this exponential growth and give the 
basic steps for deriving the Kazantsev equation, which describes the kinematic dynamo 
analytically. We derive the critical magnetic Reynolds numbers for small-scale dynamo 
action and the growth rates for different types of turbulence from the Kazantsev equation. 
In Sec. [4] we describe how the magnetic energy can be transported from the smallest scale 
of the inertial range to the largest, i. e. the nonlinear growth phase. In Sec. 5] we summarize 



our results, which are presented in much more detail in Schober et al. (2012a|b|c 



2. Basics of Our Model 



Magnetic Seed Fields. Different theories describe the origin of primordial magnetic 
fields. The first seed fields could already have been produced during inflation. , Turner fc:| 



Widrow (1988) find that Bo ~ lO^'^^-lO"^" G on a scale of 1 Mpc can be produced 



when the conformal invariance is broken. Following 'Sigl et al. (1997), there is also the 



possibility of creating a magnetic field during the phase transitions in the very early Uni- 
verse. They predict a field strength Bq w 10^^^ G from the electroweak phase transition 
and Bo ~ 10~^° G from the quantum chromodynamics (QCD) phase transition. Besides 
that magnetic fields can be generated by batteries mechanisms like the Biermann battery 
leading to typical field strengths of 10" G ( |Xu et al. |[2008l ). 

Our Approach to Turbulence. A theoretical description of turbulence starts with the 
decomposition of the velocity field v into a mean field (v) and a turbulent component 

v=(v)+(5v. (2.1) 

The correlation of two turbulent velocity components at the positions ri and r2 at the 
times t and s for a Gaussian random velocity field with zero mean, which is isotropic, 
homogeneous and (5-correlated in time, is 



{Sv,ir,,t)Sv,{r,,s)) = [(<5,, - ^) T^{r) + ^TlW <5(t - s) 



(2.2) 



with r = |ri — and the transverse and longitudinal parts of the two-point correlation 
function, Tn and Tl. Any turbulent flow can in general be described by the relation 
between the velocity v{£) and the size ^ of a velocity fluctuation. 



(2.3) 



The power-law index varies from its minimal value of = 1/3 for Kolmogorov theory 
([Kolmogor 



rov 



1941), i. e. incompressible turbulence, to Burgers turbulence (Burgers 
1948 1, i. e. highly compressible turbulence, where 1!) gets its maximal value of 1/2 ( [Schmidt 



et al. |[2009 |. 



We use the model for the correlation function of the turbulent velocity field that was 



motivated in Schober et al. ( 2012a I . The longitudinal correlation function on the different 
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length scales is 



nir) 




l) 



<r <£^ 
£^ < r < L 
L <r. 



(2.4) 



where = L Re~^/^''+^^ denotes the viscous scale of the turbulence and L and V the 
length and velocity of the largest eddies. The transverse correlation function for the 
general slope of the turbulent velocity spectrum is 
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Re(i-'')/(i+'') (£)2) < r < C 

<r < L 
L < r. 



(2.5) 



Chemical and Thermal Evolution. We determine the chemical and thermal evolution 



of gravitationally collapsing primordial gas using the one-zone model of Glover & Savin 
I ( 2009 ) . We include a modification relating the collapse time to the equation of state 
(Schleicher e t al. [2009 ) and additional Li and HeH+ chemistry ( Bovino et al. |2011b|a 



Our chemical network includes around 30 different atomic and molecular species linked by 
around 400 different chemical reactions. The initial elemental abundances are primordial 

2009 1, the initial density is no = 1 cm"^ and the 
1000 K, but we have verified that our results are 



( Cyburt 



2004 
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temperature of the gas Tq 
sensitive to these values. 

In the one-zone model the mass density p evolves as 

dp p 
dt tff 



not 



(2.6) 



where tg = ^ ii: / {?>2G p) is the free-fall time. Moreover, the temperature evolution 
determined by the energy equation, 



IS 



^^pdp_ 



de 
di 



P' 



dt 



cool 



A 



heat ; 



(2.7) 



where e is the specific internal energy, p is the thermal pressure and Acooi and Ahcat are 
the total cooling and the heating rate per unit mass, respectively. 

We find that the abundance of H is constant at low densities, but decreases at densities 
higher than about 10^° cm^'^ due to the formation of H2. For the magnetic properties of 
the primordial gas the abundances of the charged species are especially important. They 
determine for example the conductivity of the gas. At densities n < 10® cm"^, ionized 
hydrogen is the main positive ion, while at higher densities, Li+ dominates. 
Characteristic Magnetohydrodynamical Quantities. The kinematic viscosity is 



— (-) 
Ad'^n \TTm J 



1/2 



1998). Here, 



(2.8) 

p/m 



if the molecules are assumed to be rigid spheres (Choudhuri 
is the number density, d the mean particle diameter, k the Boltzmann constant, T the 
temperature and m the mean molecular weight. 

The two dominant effects that lead to dissipation of magnetic energy are Ohmic resistivity 
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VOhm and ambipolar diffusion t^ad , which we calculate according to Pinto et al. ( 2008 1 : 



„2 



?70hm,n = , (2.9a) 

47r(T||_n 



^AD,n = f ■ (2.96) 

Here cr|| crp.n and aH,n are different conductivities and the index n refers to a neutral 
species. We focus on the most important neutral species H, He and H2 and the charged 
species H"*", e~ and Li"*". The total Ohmic magnetic diffusivity is 77ohm — X^n^Ohm.n and 
the total resistivity due to ambipolar diffusion is 77AD — l/(X]n'^AD n)- 
The hydrodynamic and magnetic Reynolds numbers are defined as 

Re=— , (2.10a) 

V 

VL 

Rm= , (2.106) 

■n 

where L is the length of the largest turbulent fluctuations, i. e. the Jeans length, and V 
the typical velocity on that scale. For the calculation of the magnetic Reynolds number 
we use the sum of ryohm and 77ad ■ 

The ratio between viscosity and magnetic diffusivity is called magnetic Prandtl number 

„ Rm V 

Pm= = -. (2.11 

Re 7] ^ ' 

In Fig. [1] the density dependency of the Reynolds numbers and the magnetic Prandtl 
number is shown for both Kolmogorov and Burgers turbulence. We point out that the 
rapid decrease of Rm and Pm is caused by the dynamo amplification of the magnetic 
field. In the beginning of the collapse Ohmic resistivity is the dominant diffusion process. 
With increasing magnetic field 77AD increases proportional to and becomes the main 
process for magnetic diffusion. Since Rm and Pm are both proportional to l/rjAD, m the 
limit where 77AD 3> ??oiim, both decrease rapidly with increasing magnetic field strength. 



3. Kinematic Small-Scale Dynamo 

Phenomenology. The small-scale dynamo converts turbulent motions into magnetic 
energy. An illustrative model describing this process is the stretch-twist-fold dynamo 
( |Vainshtein fc Zeldovich 1972). The stretching of a closed magnetic flux rope leads to 
amplification of the magnetic field strength, as the magnetic flux is a conserved quantity. 
Afterwards the rope is stretched, twisted and folded such that the original shape is 
regained. This is an efficient process for increasing the field strength. 
In the limit of high magnetic Prandtl numbers this mechanism works fastest on the 
viscous scale l^, = Re"^^^''"''^^ L, as here the eddy timescale has its minimum. During 
the transition from large to small magnetic Prandtl numbers the resistive scale £jj = 
becomes larger than Then amplification takes place at roughly i^j, 
which lies within the inertial range of the turbulent velocity spectrum. Due to larger 
timescales of the turbulent eddies in the inertial range, we expect the small-scale dynamo 
to be less efficient at low magnetic Prandtl numbers. 
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Kazantsev Theory. Like the velocity field, the magnetic field can be separated into a 
mean field (B) and a fluctuation part (5B: 

B = (B)+(5B. (3.1) 

Assuming that the fluctuating component (5B is a homogeneous, isotropic Gaussian ran- 
dom field with zero mean like the velocity field, we can write down the correlation function 
as 

{5B,{vi,t)5Bj{v2,t)) = (5,, - ^) MN(r,<) + ^ML(r,i). (3.2) 
The time derivative of the correlation function is 

i ^ (fiS,) + - I ((B,) (i,,)) . (3.3) 

In the upper equation we can substitute the induction equation 

— = VxvxB-7/VxVxB, (3.4) 
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and the evolution equation of the magnetic mean field 

^ ^"^^ = V X (v) X (B) - r^off V X V X (B) 



dt 



(3.5) 



with the effective parameter rjcs = rj 
fc: Subramanian ||2005 ) this leads to 



dt 



4Kdiff 



rL(0). After a lengthy derivation (Brandenburg 



^diff 



Ml 



r 



Ml (3.6) 



with 



Kdiff(r) = ?/ + rL(o) -rL(0- (3.7) 

The prime denotes differentiation with respect to r. The diffusion of the magnetic cor- 
relations, Kdiffj contains in addition to the magnetic diffusivity rj the scale-dependent 
turbulent diffusion Tl(0) - Ti,{r). 

In order to separate the time from the spatial coordinates we use the ansatz 

ML(r,t)= ^^^(r)e2r*. (3.8) 



Substitution of this ansatz in Eq. (3.6) gives us 



d^r 



(3.9) 



This is the Kazantsev equation, which is formally similar to the quantum-mechanical 
Schrodinger equation with a "mass" h'^/{2Kdis) and the "potentiarjf] 

^ 2^ + ^ + 2(rL - Tn) 



Uir) = ''diff 



(3.10) 



2 4Kdiff 1"^ 

We solve the Kazantsev equation with the WKB approximation. For this we rewrite 
Eq. (pjl) to 



d^ 



-p{x)e{x) = 



with the substitutions ip{a 



= e^/^9{x) and 



p{x) 



U{x)]e 



2x 



Kdiff(a;) 



1 

4' 



(3.11) 



(3.12) 



where a; is a new coordinate with r = e^. Eq. (3.11) can be solved approximately by 

Vp(^da:' = |, (3.13) 

where the integration limits are the zeros of p{x). We have verified that the WKB ap- 
proximation is valid in the limit of small and large magnetic Prandtl numbers. 

Critical Magnetic Reynolds Numbers. By setting the growth rate in our equations to 
zero and solving Eq. (3.13) for Rm we calculate the critical magnetic Reynolds num- 
ber Rm^rit for small-scale dynamo action. We find that the small-scale dynamo is most 
efficient in the case of a purely rotational turbulent velocity field, i.e., for Kolmogorov 
turbulence, where Rmcrit ~ HO. The critical magnetic Reynolds number for a turbulent 
field with vanishing rotational component, i. e. , Burgers turbulence, is roughly 2700. 

f We note that there is a typo in the paper of [Schober et al~\ ( |2012a[ ), where the potential 
was derived for a general type of turbulence. The term 2Kdiff / r appeared here twice. 



lAUS 294 Dynamo & First Stars 



7 





-2000 
-4000 
-6000 
-8000 
-10000 

-12000 
20 

-20 
-40 
-60 
-80 
-100 
-120 
-140 



Kolmogorov tui'lfulence 



H \ \ [ 

Burgers turtiulence 




Figure 2. Normalized potential UL/V as 
a function of normalized distance r/L. We 
choose here Re = 10* and Rm = lO^'^ lead- 
ing to Pm = 10". 
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Figure 3. Normalized potential UL/V as 
a function of normalized distance r/L. We 
choose here Re = 10* and Rm — 10* lead- 
ing to Pm = 10"". 



In the limit of large Pm the critical magnetic Reynolds number is not really a restriction, 
as the Re needs to be larger than about 10^ for turbulence. As here Rm ^ Re, Rm 
needs to be much larger than lO'^ anyway. In the opposite limit of small Pm we have 
the condition of Rm <C Re. For low Re, Rm can fall below Rm^rit and the small-scale 
dynamo eventually can not operate. 

Growth Rates . In this paragraph we present general analytical solutions for the growth 
rate F for an arbitrary slope of the turbulent velocity spectrum i?, in the limits of large 
and low magnetic Prandtl numbers. 

In the limit of large Pm the potential has its minimum at £i, (see Fig.[2|. Thus, we expect 
the fastest growing mode to be in the viscous range, in which (3.12 1 can be approximated 
as 

P(^) = 12 (3-14) 



20 



with the definitions 



Ao = Rc(5+'5)/(2+2'') (163 - 30M) - -Re^/^T, 

3 



Bo ^ (304i9 - 98) Re^ 

,1/2 



20, 



-Re 



(3.15) 
(3.16) 



the new coordinate z = ^Re^^^Pm/3^ y and the normalized growth rate F = L/V F. 
The roots of p(z) are zi = yjE^/A^ and Z2 = ^J^^|'i R(i^'^^'^'^l'^^^+^'^ leading to a 



solution of (3.13) in the limit of large Pm 

^ (163 



- 304^) V 
60 L 



(3.17) 
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Figure 4. The growth rate of the small-scale dynamo F in the kinematic phase as a function 
of density n compared to the inverse free-fall time l/tg and the diffusion rates due to Ohmic 
dissipation and ambipolar diffusion, Fohm and Fad- 

In Fig. [5] we show the dependency of f on the Re for different types of turbulence in the 
limit of large Pm. The growth rate ranges from F cx Re^^^ for Kolmogorov turbulence 
to Re^/^ for Bur gers turbulence. Altogether we find that the growth rate increases faster 
with increasing Re when the compressibility is lower. 

The derivation of the growth rate is similar in the case of small Pm. However, the potential 
looks different in this limit (see Fig. |3| . The crucial discrepancy to the contrary limit of 
large Pm is that the potential only has a negative part between £1, and L. Thus, there 
are only real positive eigenvalues of the Kazantsev equation (3.9 1 in the viscous range. 



The resulting equations in this case are more complicated, but with the ansatz 

f = a RmT+^ 



(3.18) 



we can find an analytical solution of for the growth rate. This is motivated by our result 
in the limit of large Pm: f c>c Re''^~'^^^^^'^'^\ where we replace Re by Rm, as the latter 
determines the characteristic scale for amplification Then the pre-factor of (3.181 
turns out to be approximately 



i9(56 



103i?) .^.^ 



exp 



3 1? (56 - 103t?) 



TT (?9- 1) - 2 



(3.19) 
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with the abbreviation 



25 + \ 135i9 (56 - 103d) + (t? (79 - 157i9) ~ 25)^ - d (79 - 157?9) 
= ^ JJ^^^^^ ■ (3.20) 

We illustrate this result in Fig. [6] for different types of turbulence. 
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Figure 6. The normalised growth rate T 
as a function of Km (lower x axes) and 
Pm (upper x axes). The results shown for 
lower abscissa are only valid for small Pm, 
i. e. Rm ^ Re, while we used a fixed Re of 
10^ for the upper abscissa. 



Evolution of the Magnetic Field on Small Scales. In principle, the magnetic energy den- 



sity, Eb = i?^/(87r), evolves as 

dEB 
dt 



4 dn 
3n dt 



rohm — rAD(£^B) 



Ep 



(3.21) 



where we assume spherically symmetric collapse. By solving this equation numerically we 
find the evolution of the magnetic energy density on the viscous scale. As an initial field 
strength Bq we use 10~^° G on the viscous scale, which is a conservative value for a field 
generated by a Biermann battery ( Biermann ||1950 Xu et aZ~||2008 ). In Fig. [S] we show 
the resulting growth of the magnetic field strength. The field strength grows extremely 
rapidly as the density increases. We assume that the magnetic field is saturated as soon 
as equipartition with kinetic energy is reached. 



4. Non-Linear Small-Scale Dynamo 

During saturation on the viscous scale, the coherence length of the magnetic field shifts 



towards larger scales, a well-known behavior for the small-scale dynamo ( Schekochihin 



et al. 2002 Brandenburg & Subramanian 2005 1 , recently shown to be true also in a 



collapsing system (Sur et al. 2012). Analytical arguments suggest this occurs on the 
eddy-timescale of the current peak scale £p leading to a time evolution 



^p(t) =4(i.) 



-ft 



(4.1) 
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Figure 7. Schematic model for the time evolution of the spectrum of the magnetic field strength 
in the inertial range of turbulence. For simplicity we use a fixed frame of reference, where the 
viscous and the Jeans scale stay constant during the collapse. The different colors and line types 
represent the spectrum at different times. Further description can be found in the text. 

where ti^ is the point in time at which saturation occurs on the viscous scale. 
We illustrate our model in Fig. [7j where the different lines represent the spectrum at 
different times. The dashed green line is the spectrum at f^, the dotted red line shows 
a later time and the solid orange line represents an even later point in time at which 
the magnetic field has saturated on the Jeans scale. The slope of the curves proportional 
to is known as the Kazantsev slope in real space, which can be derived from 



the Fourier-transformed Kazantsev equation (3.9). The red line that connects the peak 
maxima at different times is a relic of the turbulence spectrum and thus is proportional 
to . Saturation on a scale I takes place, when the magnetic energy reaches equipartition 
with the turbulent kinetic energy. The saturated field strength is then 



(4.2) 



Using the Kazantsev slope, we can extrapolate the magnetic field strength onto the 
current Jeans length. By this we are able to determine the time evolution of the magnetic 
field on the Jeans scale. The result of the large-scale magnetic field is shown in Fig. |8] 
together with the field on the current peak scale and the one on the viscous scale. One 
can see that the magnetic energy is shifted rapidly onto larger scales. For Kolmogorov 
turbulence the field on the Jeans scale saturates at a density of roughly 3 cm^'^ and for 
Burgers at a density of roughly 4 cm^'^. At the end of dynamo growth we have a magnetic 
field strength of about 10~® G throughout the entire inertial range of the turbulence, 
i. e. within the Jeans volume. 

After the rapid initial dynamo amplification the only way to amplify the magnetic field on 
the Jeans scale further is gravitational compression, which leads to _B cx n?^^ for spherical 
collapse. However, the field has already reached equipartition with the kineti c en ergy at 
the end of dynamo amplification and, thus, increases only with n^/^ (see Eq. 4.2). 



5. Conclusion 

We present a semi-analytical approach to the small-scale dynamo. Our starting point 
is the Kazantsev equation which we solve within the WKB approximation in the limit 
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Figure 8. The magnetic field strength as a function of the number density on the viscous scale 
(dashed green line), the peak scale (dotted red line) and the Jeans scale (solid orange line). The 
field approaches saturation extremely rapidly, when the density has only increased from Icm"^ 
to 4cm-"^ for all types of turbulence. 



of large and small magnetic Prandtl numbers. We take into account the dependency on 
the type of turbulence, reaching from incompressible Kolmogorov turbulence to highly 
compressible Burgers turbulence. 

We determine the critical magnetic Reynolds number, which needs to be exceeded for 
small-scale dynamo action. For Kolmogorov turbulence Rmcrit ~ HO, which is in good 
agreement with previous analytical and numerical results. For Burgers turbulence we 
find that Rmcrit ~ 2700. In the kinematic phase of the small-scale dynamo the magnetic 
energy grows exponentially with a growth rate F, which appears as an eigenvalue in the 
Kazantsev equation. We find in the limit of large Pm 

rp„,»iocRe(i-'')/(i+''). (5.1) 

Here the fastest growing mode is on the viscous scale, which is determined by Re. During 
the transition to low Pm, the resistive scale, determined by Rm, becomes larger than the 
viscous one and thus the amplification takes place at larger scales. For low Pm we find 



rp^«iOcRm(i-'')/(i+''). 



(5.2) 
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Wc apply our theoretical model of the small-scale dynamo to the formation of the first 
stars. We model the chemistry and the thermodynamics during the collapse of a primor- 
dial halo and thus can determine the typical MHD quantities including Ohmic dissipation 
as well as ambipolar diffusion. We solve the energy equation and determine the evolution 
of the magnetic energy on the viscous scale. Furthermore, we model the nonlinear dy- 
namo, where we assume that the magnetic energy is shifted to larger scales on the eddy 
timcscalc. By this we are able to follow the evolution of the magnetic field on the Jeans 
scale. Diiring an increase of the density in the halo from Icm^"^ to 4cm^'^ the dynamo 
saturates quickly in terms of density for all types of turbulence. We result in dynamically 
important magnetic fields of the order of a few iiG. 
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